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Boltzmann equation for fluctuation Cooper pairs in Lawrence-Doniach model. 
Possible out-of-plane negative diff^erential conductivity 

Todor M. MishonoAQ and Yana G. ManevcQ 

Department of Theoretical Physics, Faculty of Physics, 

University of Sofia St. Kliment Ohridski, 

5 J. Bourchier Boulevard, BG-1164 Sofia, Bulgaria 

(Dated: February 5, 2008) 

The differential conductivity for the out-of-plane transport in layered cuprates is calculated for 
Lawrence-Doniach model in the framework of time- dependent Ginzburg-Landau (TDGL) theory. 
The TDGL equation for the superconducting order parameter is solved in the presence of Langevin 
external noise, describing the birth of fluctuation Cooper pairs. The TDGL correlator of the super- 
conducting order parameter is calculated in momentum representation and it is shown that the so 
defined number of particles obeys the Boltzmann equation. The fiuctuation conductivity is given by 
an integral over the Josephson phase 6 of the particles distribution, depending on that phase, n(6) 
and their velocity v{9). It is demonstrated that in case of overcooling under Tc the transition from 
normal to superconducting phase while reducing the external electric field runs via annulation of 
the differential conductivity. The presented results can be used for analysis and experimental data 
processing of measurements of differential conductivity and fluctuation current in strong electric 
and magnetic fields. The possible usage of a negative differential conductivity (NDC) for generation 
of THz oscillations is shortly discussed. 

PACS numbers: 74.40.-|-k, 74.20.De, 74.25.Fy, 74.72.-h 



I. INTRODUCTION 



Fluctuation phenomena in high-temperature cuprates 
(CUO2) are greatly pronounced because of the short 
length of coherence. The number of degrees of freedom in 
the sample's volume is V/^^^(0)i^c(0). The fluctuations in 
strongly anisotropic superconductors are approximately 
two-dimensional in layers at a distance s from each other 
and the number of fluctuation modes is V/^Qf,(0)s. A re- 
duction of the dimension, as a rule, leads to enhance- 
ment of the fluctuation effects. The paraconductivity of 
Aslamazov-Larkin (AL), when to the normal conductiv- 
ity we have to add the contribution of the fluctuation 
Cooper pairs is the most thoroughly investigated fluc- 
tuation phenomenon. The AL conductivity is easier to 
be observed if the normal phase conductivity is small. 
That is the case of a strongly disordered conventional su- 
perconductor or high-temperature cuprates, containing 
CUO2 planes as a main structural detail. The amplitudes 
of the electron hopping between Cu02 layers are small, 
which leads to intensive diminution of the conductivity 
in the so called dielectric z-direction, perpendicular to 
the Cu02 layers. As a consequence one could apply con- 
siderable electric fields that will not result in significant 
heating of the sample. The small heating gives the oppor- 
tunity to investigate how the strong electric field's nonlin- 
ear effects influence the fluctuation conductivity a. Such 
investigations could reveal new details in the kinetics of 
the superconducting order parameter. They could spec- 
ify the parameters of time-dependent Ginzburg-Landau 
(TDGL) theory which describes, with acceptable accu- 
racy, the fluctuation phenomena in almost all supercon- 
ductors. We recommend the monograph by Larkin and 
Varlamov- for a general review on fluctuation phenomena 



in superconductors and the review article (and references 
therein)^ which is especially devoted to Gaussian fluctua- 
tions in layered superconductors. In addition, the recent 
experimental investigations on cut-off effects^ and reduc- 
tion of paraconductivity with increase of the electric field 
in cuprates^ are also recommended. 

The purpose of the current work is to demonstrate that 
negative differential conductivity (NDC) can be reached 
in the specific nonequilibrium conditions of a supercooled 
below Tc superconductor in normal state in external elec- 
tric field Ec- This NDC can find possible applications in 
the usage of layered cuprates as active media in THz 
generators. We derived explicit formulas for differential 
conductivity in out of plane direction and for the current 
in case of external magnetic field B^ ■ In the special case of 
evanescent magnetic field our results agree with the work 
by Puica and LangJ^ Our TDGL result for the momen- 
tum distribution of the fluctuation Cooper pairs obey the 
Boltzmann equation. The new results for the current j^ 
can be used for the experimental data processing and in- 
terpretation of the experimental results for out-of-plane 
fluctuation conductivity for samples in the geometry of 
plane capacitor with a layered superconductor between 
the plates. Special attention is paid to the analysis of for- 
mulas for the differential conductivity in different phys- 
ical conditions. How to reach easily the regime of high 
frequency oscillations is briefly considered. 

For this sake the time dependence of the order param- 
eter is estimated under Langevin approach and its corre- 
lator in coinciding time arguments is calculated. We have 
demonstrated that this correlator in case of homogeneous 
electric field obeys the master equation. In this way 
Boltzmann equation for the fluctuation Cooper pairs has 
been derived, which solution gives the momentum dis- 



tribution of the superconducting order parameter. The 
electric current is represented as an integral over this dis- 
tribution and the velocity. The derived formula for the 
current could be generalized for a self-consistent treat- 
ment of the interaction between the fluctuations. Such a 
possibility has been discussed in short. We have shown 
that when T < Tc the phase transition to the supercon- 
ducting state passes through the critical point of zero 
differential conductivity crdiff' = 0. This qualitative result 
is one of the important consequences predicted by the 
theory which can be described by the deduced formula 
for the current jz(e, V , h) as a function of the dimension- 
less electric field V — eEzTos/h, dimensionless magnetic 
field h = 27r^^j(0)i3z/$o and the reduced temperature 
e = ln(r/rc) « (T - Tc)/Tc Eq. §5\). Further decrease 
of the electric field E in the regime T < Tc could originate 
either electric oscillations or lead to formation of domain 
structure, but for sure the homogeneous phase loses sta- 
bility and a new physics takes place. What we do is to 
offer systematical measurements of the generation of har- 
monics for investigation the role of strong electric fields. 
The theoretical analysis is completed in the framework 
of the TDGL theory applied to the Lawrence-Doniach 
(LD) model for a layered superconductor, introduced in 
the next section. 



II. MODEL 

Our starting point is the TDGL equation, which in 
momentum space takes the form 

[e{P-e*A)+aoe]iPp{t)=-2aQTo[dtijp{t)-Cp{t)], (1) 

where the argument of the kinetic energy is the gauge 
invariant kinetic momentum 



p{t)=P~e*A{t), 



(2) 



t is the time, and tq is the TDGL relaxation time. The 
canonic one P — const is a conserving quantity for a free 
particle moving in an external electromagnetic field. In 
LD model the energy spectrum is given by 



£(P) 



pI+pI 

2mab 



aar 



(1-008 6*). 



(3) 



Here we have used standard notations for the time- 
dependent vector potential of the homogeneous electric 
field perpendicular to the layers in our analysis Ez(t) = 
—dtAz{t), the charge of the Cooper pairs |e*| = 2|e|, their 
effective mass in the plane of the layers for the uniaxial 
superconductors which we study iriab, the effective mass 
for motion in c-direction rric and the GL parameter 



Oo 



aor 



2mabeab{0) 2m,e?(0)' 



nicS^ 



(4) 



For high values of the LD parameter r = (2^c(0)/s) 
their model simply describes the anisotropic GL theory, 



applied to an uniaxial crystal {nia = rUb ^ rric)- Con- 
versely, when the coherence length in perpendicular di- 
rection ^c(O) is much smaller than the distance between 
the layers we have a system of two-dimensional supercon- 
ductors with a Josephson coupling. 
The gauge invariant Josephson's phase 



eit) . ^ 



Pz-e*Az{t) 
h/s 



e (-TT, +7r) 



(5) 



actually is the kinetic momentum of the Cooper pair mea- 
sured in units, connected to the lattice constant in z di- 
rection s. For the sake of simplicity from now on we 
are going to omit the z index, which would be under- 
stood by the context. For the relaxation time tq in the 
TDGL equation, the Bardeen-Cooper-Schrieffer (BCS) 
theory gives 



H 

2^ 



■n 



(6) 



This result for the microscopic theory is extremely ro- 
bust. The same value we have for dirty bulk supercon- 
ductors with isotropic gap and clean two-dimensional d- 
type superconductors i Experimental researches in high- 
Tc cuprates confirmed this value within the experimental 
accuracy. For a discussion of different methods for de- 
termination of the relaxation time tq, and references of 
relevant experimental works see Sec. 4.2 and Sec. 4.3 
of the review Ref. ^. That is why TDGL theory does 
not need modifications when applied for d-wave layered 
cuprates. 

The GL theory is applicable for small by absolute value 
reduced temperature |e| ^ 1, whereas we are going to use 
the negative values of this parameter for a description 
of an overcooled superconductor plugged in an external 
electric field. For the Langevin term in the TDGL theory 
we have a white noise correlator 



T 



(Cpi(^l)Cp2(i2)) = 5p,^pJ{ti-t2) 



aoTo 



(7) 



Analyzing bulk properties of the material we formally 
assume periodic boundary conditions for a great distance 
i ^ s, so that the quasi- momentum takes the discrete 
form 



27r?i / T:h T:h 
Pz = ^—(integer) G , H 

L \ 3 s 



(8) 



Let us introduce convenient for this task dimensionless 
variables: u = t/r for the dimensionless time and 

^(k, e)^'-^ = kl + kl + Ul- cosO) (9) 

ao 2 

for the dimensionless kinetic energy, where 

, _ P.u^) , _ Pvib{Q) , ^ PzUo) 



2fc./v^, 



(10) 



are the components of the dimensionless wave-vector 
and the Josephson's phase. For the dimensionless noise 
(*{u) = To(p{t) the correlator looks like 



(C(wi)Cp("2)) = n^S{ui - U2), 



T 



(11) 



For the sake of simplicity let us first consider the 
one-dimensional case of hopping of fluctuation current 
through a chain of Josephson junctions. Then the TDGL 
equation in terms of dimensionless variables reads 



du'ipqiu) = -iy{u)'lljq{u) + Cq{u), 



(12) 



where q = Pz^c{^)/fi is the conserving dimension- 
less canonical momentum in z direction, A{u) = 
— e* A{u)£^c{^) / h is the dimensionless potential momen- 
tum corresponding to the dimensionless electric field 

f{u) = e*E{t)ToiMlh = d^A{u). (13) 

The dimensionless decay rate is given by 

2v{u) = uj{q + A{u)) + t (14) 

and the kinetic dimensionless momentum for the one- 
dimensional task we are solving is 



kz=q + A{u). 



(15) 



III. DERIVATION AND SOLUTION OF THE 
BOLTZMANN EQUATION 



The TDGL equation Eq. (|12p is a linear ordinary differ- 
ential equation and one can easily check that its solution 
has the form 



V'g(M) 



C,q{u2) exp 



«2=0 

■0,(0)} exp 



v{ui)(lui 

ui=0 

v{u-i)du-i 



U2 



«3=0 



(16) 



The averaging of the superconductor's order parameter 
over the noise gives the distribution of the mean number 
of particles with respect to the canonical momentum 

nq{u) == (V'*(u)V'g(M)) 



n^ I du2 exp 
'0 

+ nq{u = 0) exp 



'2 I dui i/(Mi) 

'U2 



v{us) dw3 



(17) 



In case of constant electric field the vector-potential and 
the decay rate take the form 

A(7i) = /u, 2v{ii) = - [1 - cos((/) + 2Vu)] + e, (18) 



where we have introduced the notations 



/ eE^TQS 



PzS 



^r 



h 



Tr 



const. (19) 



An elementary integration of the upper equation (|17p in 
the limit u — > 00, using the trigonometric relation 



sin((/) + 27'u) — sin = 2 cos (0 + Vii) sin {Vu) 



gives 



n{e) 



du 



X exp 



^+2.'" 



— sin(7'u) cos [9 - Vu) . (20) 

Here we have turned to a new integration variable 
u = u — U2, see Eq. p7)) . which describes the birth of 
fluctuation Cooper pairs in the moment tqu before the 
time of observation t ~ tqu. In such a case the integrant 
gives the age distribution of fluctuation Cooper pairs. 
The periodical dependence sin(0 + 2Vu) on time m is a 
reflection of the Bloch oscillations of the free charged par- 
ticles which move in a periodical potential and constant 
electric field. 

In case of weak electric fields s\\\{Vu) « Vu we obtain 
the equilibrium distribution density function 



n{e)^ 



(1 -COS 6*) -he' 



(21) 



The reconstruction of the plane components of the mo- 
mentum leads to the general formula of Rayleigh-Jeans 
for the equilibrium distribution of the particles 



T 



T 



2m„(, 



+ _^ (1 _ cos %^) 



floe 



e{p)- n 



(22) 

Hereby formally we may identify the chemical potential 
with the Landau parameter fi = —aoe. In this sense the 
critical temperature Tc sets the point of annulation of 
the chemical potential and fixes the beginning of Bose- 
Einstein condensation of the fluctuation Cooper pairs. It 
may be easily checked that the momentum distribution 
of the particles obeys the master equation of Boltzmann 



dtnp(i) 



np(i) - np-e'Ait) 



Tp-e'A 



npjt) 

TP-e-A To 



(23) 

where the decay rate is proportional to the kinetic energy 
measured from the chemical potential 

^— = 1. (P - e*A(i)) = £M±^. (24) 

Tp-e'A To ao^o 

In such a way we derived the Boltzmann equation for 
fluctuation Coper pairs for out-of-plane transport of a 
layered superconductor. The argument of the distribu- 
tion and decay rate is the time-dependent kinetic mo- 
mentum. The direct derivation of Boltzmann equation 
directly from the TDGL theory is too long and contains 
a lot of technical details.— 

In case of zero electric field the relaxation rate of the 
order parameter with zero value of the wave vector runs 
through a typical slowing down 



T{e) 



~p = 



e 



(25) 



This is the characteristic time of "drying" at e > of 
the spatiaUy homogeneous Bose condensate with an or- 
der parameter ijj{r) = const. If we consider the distri- 
bution depending on the kinetic momentum we obtain 
from TDGL theory the standard form of the Boltzmann 
equation derived in 1872 

dtn{p,t) + e*E{t)-dpn{p,t) 

n{p,t)-nip) n{p,t) n^ 

t(p) t(p) To 

In our dimensionless variables for the static case this 
equation reads 



2Vden{0) 



-(l-cos6')-|-e n{e) + 



(27) 



which solution is Eq. ([20|) . 

In such a way we derived the Cooper pair Stoss-integral 
in the framework of the TDGL theory, the result is for- 
mally an energy dependent r- approximation with a con- 
stant birth rate in momentum space 






To 



(28) 



This integral describes collisions between normal parti- 
cles and creation of Cooper pairs and back process of 
the decay of the condensate particles. We wish to point 
out that the TDGL equation is a diffusion-like equa- 
tion which does not lead to quasi-classical dynamics of 
quasi-particles.^ This shows that Boltzmann equation 
has broader applicability than arguments used in his text- 
book presentations. 

Let us now turn back to the analysis of the distri- 
bution function. In the continual limit when the co- 
herence length significantly exceeds the lattice constant 
^c(O) 3> s or in other words when the LD parameter is 
big enough r ^ 1 we have to consider small Josephson's 
phase \6\ ^ 1 and weak fields {Vul <C 1. The distribution 
over the kinetic momentum Eq. (|20p then becomes2.i§ 



n(k^) 



exp 



^{e + kl)u + kju^ 



1 



f^u-" ^ du. 

(29) 

While in Refs. 130 are analyzed bulk superconductors in 
the present work we consider a layered superconductor 
with Josephson coupling between the layers for which we 
apply LD model. The last formula may be derived di- 
rectly from the Boltzmann equation applied in case of 
parabolic dispersion e{pz) ~ p1/2mc. After we have in- 
vestigated the momentum distribution here in the next 
section we will obtain a general formula for the density 
of the electrical current. 



IV. FLUCTUATION CURRENT 

Let us start our analysis with the one-dimensional case 
as the results may be easily extrapolated for the case of 



higher dimensions. The density of the current 

Pz 



J; 



(30) 



is a product of the particles' charge e*, their density 
n{pz)/ L and velocity 



de fi . , 

v{Pz) = ^— = sm^ 

opz rricS 



(31) 



This formula is a special case of the general procedure 
for derivation of the current density in the framework of 
the GL's theory (Ref. @, sec. 115; Ref. [13], sec. 45) 



J(r) 



6A{r) 



^*(r)e(-iW-e*A(r))V'(r)d^x, (32) 



where after the functional integration the vector poten- 
tial is put to be spatially homogeneous and the order 
parameter takes the form of a plane wave 



|2j3 



V'(r) = J^e'P-"^'', np= |V(r)|^d^x, V = / d^x. 

(33) 
In the one-dimensional case, imagine current along a 
chain of Josephson junctions, if we use the relation 



E 



2Trh 



dpz 



de 

2^' 



(34) 



the substitution of the Boltzmann equation's solution 
Eq. ([20| in the general formula for the current Eq. ((30|) 
gives 



,-(iD) - 



e*rT 



sin Me* 



■{-{^ 



cxp <; - ( e + -J u 



— — sm{Vu) cos {9 — Vu) f du. 



i)}: 



(35) 



Here we have taken into account that 9 — Pz-s/h and 
aon^ — T, see Eq. ([5]) and Eq. pT|) . The averaging with 
respect to the Josephson's phase 9 can be expressed by 
the modified Bessel functions^^ 



/ c\9 
— cos m9e 



Thus the formula for the current reads 



,-(iD) 



e*rT 
2h 



due-(^+'~/2)fi s[n{ru)Ii 



(36) 

rsai{Vu) 
2P 



(37) 
In order to calculate the out-of-plane component of the 
current jz in the LD model we have to consider the com- 
plementary decay rate excited by the kinetic energy in 
the xy plane, i.e a6 Cu02 plane. For this reason formula 
Eq. ([M]) has to be modified as follows 



2.= £M±^!l(!)^!:(l_cos0)+, + e, (38) 

oo 2 



where according to Eq. ([TO 



Px + Py _ . 2 , . 2 

Ztriabao 



(39) 



is the dimensionless in-plane kinetic energy Hmited by 
a cut-off parameter c; for more details see, for example, 
Ref. I^. We have to sum over the kinetic momentums in 
the a6-plane 



1 



dpxdpy J pI+pI \ ^ 
{2TThf\2mabaoJ 4^eab(0) Jo 



dvfiv)- (40) 



This additional integration with respect to the two- 
dimensional degrees of freedom taken into account in 
Eq. ((571) gives the replacement 



^<biO) Jo 



-rju 



Arj 



e-<^"(l-e™) 



(41) 



Then the expression for the fluctuation current takes the 
final form 



J.(e"^) 



erT 



4^^e^fc(o) Ju^o 



dSe-(^"+'^/2)« 



1-e- 



X sin(7'u)/i 



rsin(7''u) 



(42) 



where e means the self-consistent treatment of e param- 
eter, which will be analyzed in the next section. 



V. MAXWELL-HARTREE SELF-CONSISTENT 
APPROXIMATION 

Up to this moment we have neglected the influence 
of the nonlinear term in the GL theory. The coefficient 
standing in front of this term participates, for instance, 
in the equation for the equilibrium value of the order 
parameter above Tc 



[aoe-f^lV-l'lV-^O. 



(43) 



The idea of the self-consistent approximation (SCA) is 
to replace the density of particles with its value aver- 
aged over all fluctuation modes. In such a way we obtain 
a self-consistent equation for the renormalized reduced 
temperature, see Ref. |j. Sec. 3.3 



= In ■ 



T b 



,i3D) 



(e-)- 



(44) 



The idea is coming from the Maxwell treatment of the 
ring of Saturn - the first work on collective phenomena in 
physics. Maxwell concluded in 1856 that the ring cannot 
be a rigid object, but consists of "indefinite number of 
unconnected particles" and each of them is moving in the 
averaged gravitational field of the others. Analogously 
the fluctuation density of Cooper pairs renormalizes its 
chemical potential as shown in Eq. (|44p . 



At flrst from Eq. ((20)) we calculate the one dimensional 
density of the Cooper pairs 



,(1D) = f ^„(0) 



(45) 



Analogously to the calculation of the current Eqs. ((5D|) 
and (1341) now we have 



s Jo 



2V 



The volume density can be obtained via summation of 
this one-dimensional density with respect to the in-plane 
degrees of freedom 



n(3D)(e-) 



dpxdpy (ID) /- pI+pI 



{2TrhY 



2ml^ao 



(47) 



cf. the analogous transition between 2D and 3D case for 
in-plane conductivity of layered superconductors, Ref. y, 
Sec. 2.3. The kinetic energy cut-off, Ref. 0, Eq. (11) 



pI+pI 

2m* 



< aoc 



(48) 



leads to the same replacement as in Eq. ((4T|) and for the 
3D density we get 



n^'°He,r) 



ut 



,1-e- 



4Tsea6(0) A=o 
xe-(^"+'-/2)-/„ ( 



rsinP 
2V 



(49) 



Additionally we may represent the nonlinear parame- 
ter h in terms of the Ginzburg-Landau parameter kql = 
Aa6(0)/Cab(0) and the flux quantum $o, see Ref. 0, 
Eq. (165) 



b-2Jp^]' 
K'^omlJ 



(50) 



We also introduce the dimensionless Ginzburg number 
for a layered superconductor (cf. Ref. Q, Eq. (168)) 



eGi(T) = 2Trfio — 



T f Xab{0)Y _ 2fio^l^e^eaMT 



s V $0 



nfi^s 



(51) 



In terms of the so introduced notations the self- 
consistent equation for the renormalized temperature 
reads 



e = e + CGi / dw 



1-e- 



^"[ 2V J' 
(52 



in agreement with the results of Puica and Lang Ref|5|, 
Eqs. (11) and (17), where the temperature-dependent 
Ginzburg number eGi(T) is denoted by gT. The solution 
e of equation ((52p has to be substituted as an argument 



in the formula for the current Eq. ([42)1 . This agreement 
demonstrates that it is possible to derive in the begin- 
ning an 1-dimensional (ID) formula for the Josephson 
chain and later on to generalize this approach to a 3- 
dimensional (3D) case of layered superconductor. In the 
next section we will analyze the modification of this re- 
sult in case of out-of-plane external magnetic field. 



VI. INFLUENCE OF A PERPENDICULAR 
MAGNETIC FIELD. MAGNETOCONDUCTIVITY 

A. General case 



In particular, in the exponential dependence, which we 
have in Eq. (|4ip , we have to sum up the limited geometric 
progression 



/" 

Jo 



N^-l 



d?7e-'"~'~>2/i Y^ e" 



■{2n+l)hu 



(59) 



n=0 



which more easily can be represented as a subtraction of 
two infinite ones 



Na-l 



OO CXD 



E-E-E 



(60) 



n=0 n=0 



.N, 



Let us consider an external magnetic field also applied 
perpendicularly to the Cu02 planes. In this case the vari- 
ables can be separated and the task for the calculation 
of the current reduces to the ID case Eq. p7|) . which 
we have already considered. For the two-dimensional 
movement the magnetic field arouses equidistant oscil- 
lator spectrum 



Pab 

2m* 



huJr 



e*B, 



(53) 



The degeneration rate is determined by the Landau sub- 
zone capacity -B/$o, where $0 = 2TTh/e* is the magnetic 
flux quantum. For parametrization of the GL theory it 
is convenient to introduce the linear extrapolated upper 
critical field, represented via the coherence length in the 
a6-plane 



Bc2i0) 



-Tr 



Bc2{T) 



dT 



$0 



Tc 



2<,(0)- 



(54) 



For facilitation we will work with the dimensionless mag- 
netic field 



B, 



huJr 



Bc2{0) 2ao' 



(55) 



with the help of which the dimensionless spectrum quan- 
tizes as 



■q^ {2n+l)h, 71 = 0, 1, 2, 3, 



(56) 



The maximal kinetic energy is reached at some big num- 
ber Nc, which cuts off the summation over the Landau 
levels 



hujcNc = OqC, 



2hNr 



(57) 



The influence of the magnetic fleld reduces to the re- 
placement of the integration over the kinetic energy to a 
restricted summation over the Landau levels 



No-l 



r drjfirj) ^ Ar^ Y. f ii^ri + l)h) 
■'0 „=o 



(58) 



In this way the summation over the in-plane degrees of 
freedom simply reduces to the shift 



dTye""'^ 



sinh hu 



1-e- 



(61) 



Then the substitution described in Eq. (|¥l]) in the pres- 
ence of external magnetic field takes the form 



e--d,^SfeS^. (62) 



47rCifc(0) Jo 



4<ah(0)sinh/iw' 



This is the recipe with which one may switch from the 
already solved one-dimensional task Eq. p7|) to the three- 
dimensional. Comparison with the case of zero magnetic 
field shows that we may use the 3D formulas Eqs. ([5^ 
and (j42p . where we only have to make the substitution 



h 



u smh{hu) 



(63) 



In this way, in agreement with Ref. 5, Eqs. (13) and 
(14), we derive the final formulas for the self-consistent 
reduced temperature 



e == e + ec 



Jq smh hu \ 



am{Vu) 



2V 



(64) 

and the perpendicular to the Cu02 planes density of the 
current 



Jz{e,V,h) 



16fts^afc(0) A=o sinh/iu 

s\n{Vu) /rsiniVu 



-h 



V \ 2V 

where we have used the relation 

erTV e^il(!S)E, 



^^KaM 16/i<^,(0) 



(65) 



(66) 



This explicit formula is one of the new results of the 
present work. Let us now analyze the magnetoconduc- 
tivity in some characteristic particular examples. 



B. Weak magnetic fields 

Let us consider the application of the upper formulas 
in case of weak magnetic field /i <C | e| . If we use the well- 
known summation of Euler-MacLaurin (cf. with Ref. 0, 
Eq. (23)) we have 
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For weak magnetic fields we may use the approximation 
h 11 



sinh hit u 6 



-h^u. 



(68) 



The substitution of this decay into Eq. (|65p leads to the 
special expression for the current in the presence of weak 
magnetic fields 
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The coefficient in front of the integral over the reducing 
term may be expressed via the coherent lengths and the 
lattice constant as follows 
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The experimental data fitting with this formula could 
lead to more precise specification of i^c(O) for the investi- 
gated sample. The procedure Eq. ([68|) applied to Eq. (jM]) 
ends in a self-consistent equation for the reduced temper- 
ature in case of weak magnetic fields 
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In the next subsection we will analyze the case of values 
of the magnetic field close to the phase curve Bc2{T), 
where the influence of the fluctuation pairs is most vividly 
expressed. 



In the case of strong magnetic fields it is more appropriate 
to use the ID formula Eq. (|57|) . where we have to add a 
summation over the Landau levels, taking into account 
the degeneration and the cut off 

oo 

, (3D) (- V,h)^^Y^ j^^^ (e. + 2nh, V). (73) 

" n=0 

In close vicinity of the upper critical fleld |e| <C ^ in the 
summation the share of the lowest Landau level domi- 
nates over the others and its contribution gives 
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(74) 

Similar procedure applied to the number of particles' 
density Eq. ([46]), see also Eq. ((52|l . 



e/i = e-h 



bB 



n^''\eH,V) 



(75) 



determines the equation for the self-consistent reduced 
temperature 

e. = .. + 2/..ci^°"d^e-(^-^-/^)-/o ((-^^) • (76) 

Close to the upper critical magnetic fleld the fluctuations 
are practically one-dimensional and such a reduction of 
the space dimensionality signiflcantly amplifies the fluc- 
tuation phenomena. 

Next we are going to derive an explicit formula for the 
differential conductivity in the most general case of LD 
model and analyze it in the particular cases of both weak 
and strong magnetic flelds. The physical conditions of 
possible occurrence of negative differential conductivity 
will be discussed. 



VII. DIFFERENTIAL CONDUCTIVITY. 
POSSIBLE NDC 

The non- linear effects over the conductivity azz{Ez) = 
jz/Ez could be best understood in the investigation of 
the differential conductivity 






(77) 



C. Strong fields close to B^2{T) 

Close to the phase curve it is more convenient to ac- 
count the reduced temperature in regards of the phase 
transition temperature, which is a function of the ex- 
ternal magnetic fleld Tc2{B). Comparison with Eq. ([M|) 
reads 
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(72) 



In order to derive it we have to differentiate the gen- 
eral formula for the current in an external magnetic field 
Eq. (j65p and use the well-known properties of the modi- 
fied Bessel functions 



^ dl„i{z) 
dz 



l{z) + Ir,i-l{z) 



(78) 



which can be found in any special functions manual, for 
example see Ref. [Ul or Ref. |T2 and references therein. We 
turn to differentiation with respect to the dimensionless 
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electric field P, taking into account the value of the decay 
rate tq that follows from the BCS theory 

dV esTn |-Rr;s^ tt/i 
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(79) 



From Eq. (|65p it is straight forward that the electric field- 
dependent part in the differential conductivity is given by 
the expression 
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where for the sake of brevity we have introduced the no- 
tations 



s — sm(Vu), c = cos(7'u), /,„ = /„ 
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(81) 
Having in mind the considerations above it can be easily 
shown that the general formula for the differential con- 
ductivity takes the final form 
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This is the equation we are going to use in order to ana- 
lyze the presence of negative differential conductivity for 
given samples. Hereof we may derive the expressions for 
(T(jiff in the special cases of weak and strong magnetic 
fields. 

In the absence of external magnetic field the limit h —> 
gives h/ sinh h —^ 1 and we obtain 
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which may be derived straight from Eq. 

For strong magnetic fields we take into account only 
the contribution of the lowest Landau level and the dif- 
ferentiation of Eq. (|7i|) reads 
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If we want to return to the one-dimensional differen- 
tial conductivity for strong magnetic fields close to the 
phase curve we only have to omit the number of magnetic 
threads B^/^o in the upper equation, which gives 
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Of course the calculated derivative {djz/dEz)e does 
not describe the whole effect. We have to take into ac- 
count also {de/dEz)e and the increase of the sample's 
temperature above the ambient temperature 
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due to the thermal resistance. The "chemical" reaction 
of pairing of normal charge carriers e + e —> e* creates 
also a decreasing of the number of normal charge carriers, 
and in the Drude formula, for example, we have to take 
into accounti the density of state's corrections (DOS) 

CTnorm = ("e " 2n{e,V))e'^Tnorm/me- 

Al those effects are however smaller compared to the 
dominating Aslamazov-Larkin conductivity which quali- 
tatively can describe the appearance of new physical ef- 
fects such as NDC, for example. 

In order to better understand the nature of NDC in 
supercooling regime in the next section we will analyze 
the Green functions of the Boltzmann equation. 



VIII. GREEN FUNCTIONS OF THE 
BOLTZMANN EQUATION 

Now let us paint some Cooper pairs in green in order to 
trace their motion i.e. in order to analyze the influence of 
the random noise on the momentum distribution of the 
Cooper pairs we will replace the constant income term in 
the Stoss- integral hq/tq with a (5- function inhomogeneous 
term in Eqs. §^ and ^ 
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For simplicity we will consider ID case. 

Before the moment of creation there are no Cooper 
pairs np{t < 0) = 0. Then immediately after the influence 
we have ^-like initial distribution 

2Trh. 



np{t = +0) = Sp 



-5{p-po). 



This distribution is normalized to have one Cooper pair 
in the initial moment 
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born by thermal fluctuations; pictorially speaking "by 
the sea foam." According to the separation of variables 
in the Boltzmann equation Eq. (|23p in every momentum 
point the evolution of the distribution function is inde- 
pendent 
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The substitution here of the decay rate from Eq. (fTSl) 
gives 
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It is instructive to consider the 3D limit case of |e| ^ r; 
formally one can consider the unrealistic for the cuprates 
case of r 3> 1 . For small angles 
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the Taylor expansion of the trigonometric functions in 
Eq. dH]) gives 
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The simplest example is to consider zero initial momen- 
tum q = S,c{0)po/h = 0. In the supercooling regime where 
e < in the beginning the number of Cooper pairs in- 
creases. It is a typical lasing process - increasing of the 
number of coherent Bose particles in one mode; imagine 
a drop of rain increasing by the condensation in a humid 
atmosphere. Finally, however, the ^f'^u^ term dominates 
in the argument of the exponent and Af{oo) = 0; one can 
say that the droplet is evaporated by the big kinetic en- 
ergy like a meteorite falling in the earth atmosphere. The 
droplets are not smeared during their lifetime 
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They have just a drift in the kinetic momentum space. 
This drift related to the electric current is analogous to 
the rain fall created by the earth acceleration. 
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In upper Green function solution of the Boltzmann equa- 
tion Eq. ((26)) we have recovered the 3D variables. 

Now we can understand qualitatively the mechanism of 
creation of NDC for small electric fields in a supercooled 
superconductor. Let us take for illustration po = case 
in Eq. ([93l). For e < 0, ft- = 0, and |e| < r this function 
has a maximum^ 
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The maximal increase of the "weight" of the Bose droplet 
A/'(Minax)/A/'(0) can reach big values for small enough 
electric fields 
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Decreasing of the electric field increases the maximal 
size of the Cooper pair droplets and the electric current 
which is proportional to the number of particles in the 
droplets. This precursor of the Bose condensation cre- 
ates the NDC. Our preliminary numerical calculations 
for YBa2Cu307-5 (Y:123) taking the parameters from 
Ref. iJ: s = 1.17 nm, Cafc(O) = 1.2 nm, ^^(0) = 0.14 nm, 
Te = 92.01 K, ri"™^ = 86.9 K, kgl = 70, l/cri"°™^ = 
5 milcm, c = 0.5 demonstrated existence of NDC; this 
gives the hope that NDC can be experimentally observed 
for out-of-plane transport in Y:123. 



IX. ANALYSIS, DISCUSSION AND 
CONCLUSION 

Let us start with the formal discussion. In case of big 
enough overcooling e + r/2 < 0, when e < and |e| > r 
the integral for the current Eq. (|42|) is divergent. This 
means that the electric field cannot prevent the appear- 
ance of a superconducting condensation, as is the case of 
volume superconductors. The reason for this is that the 
one-dimensional zone in the energy spectrum Eq. ([3]) has 
a finite width aor, in contrast to the LD model, which 
shall be obtained if we take the continual limit of r ^ 1. 
However, if the LD parameter is small enough we have to 
take into account the interaction between the fluctuation 
modes, described by the nonlinear term in the GL the- 
ory. These effects become important when the reduced 
temperature e is comparable to the Ginzburg number for 
layered superconductors 



eGi(Te) = 



1 



4<\(0)sAC' 



(100) 



where (see Ref.0, Sec. 3.3, Eq. 169) AC is the jump of the 
heat capacity, extrapolated from the critical behaviour of 
the heat capacity C{T) as a fitting parameter, or can be 
evaluated by the electrodynamic properties 
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In the Maxwell-Hartree SCA the total volume density of 
the fluctuation Cooper pairs leads to an effective heating 
(see Ref.i, Sec. 3.3, Eqs. (164) and (173)). As a whole, 
inclusion of the SCA does not change the qualitative be- 
havior of the obtained results. In overcooled supercon- 
ductors e < 0, for example, we expect the appearance of 
a negative differential conductivity (NDC) with the de- 
creasing of the constant electric field. Whether this will 
lead to formation of a domain structure of normal and 
superconducting layers or to existence of electric oscilla- 
tions is a question of further analysis. Our theory, how- 
ever, predicts that this phase transition when T < Tc 
will take place via annulation of the differential conduc- 
tivity. In this respect the layered superconductors give 
an advantage as, because of the strong anisotropy, the 
normal conductivity perpendicular to the layers is too 
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small. Thus the sample will not heat intensively, which 
gives the opportunity to be observed the nonlinear effects 
of a strong electric field over the fluctuation conductivity. 
These nonlinear effects can be observed as an amplifica- 
tion of the current harmonics generation as we reach the 
critical region. Systematical investigation of these har- 
monics will lead to specification of the parameters in the 
TDGL theory and to clarification of the kinetics of the 
superconducting order parameter. Here we would like to 
consider the possible application of our theory for dif- 
ferent layered cuprates. For instance, for the extremely 
anisotropic Bi2Sr2CaCu208-i-£c (Bi:2212) the value of the 
LD parameter r is exceptionally small. The heating from 
the electric field applied in the "dielectric" z-direction is 
very weak, but the area where current harmonics can be 
observed is very narrow as well. Because of this their ob- 
servation would require samples of extraordinary quan- 
tity. Wider and easier for observation would be the non- 
linear region for the cuprates with a moderate anisotropy 
like Y:123. For most favorable we consider the cuprate 
Ta2Sr2CuOe (Ta:2201) where a coherent Fermi surface is 
observed. When the anisotropy is small and r parameter 
bigger, the half-width of the z-zone |aor = h^ /rn^s^ is 
bigger and then the electric field E^ causes more signif- 
icant increase of the kinetic energy of the Cooper pairs 
200^(1 — cospzs/h). Because of this, for instance, the 
overcooling region when e < and V ^ would be wider 
and more easily observed. Exactly, in such more easily 
reached by the experiment conditions, wc hope to observe 
the annulation of the total differential conductivity. For 
this purpose submicron high qualitative films of Ta:2201 
are necessary to be used. The investigation of the har- 
monics generation above T^, is a useful beginning for the 
systematical research of the nonlinear effects of the strong 
electric field over the fiuctuation current in z-direction. 

The experiment which we advocate is principally very 
simple, but according to the best of our knowledge not 
done yet. The fluctuation conductivity have to be investi- 
gated in a voltage biased circuit; a constant (DC) voltage 
has to be applied. A significantly smaller AC voltage has 
also to be applied to the sample. The DC current gives 
the conductivity and a measurement with a smaller AC 
current gives the differential conductivity. The circuit 
should be similar to the devices invented for investigation 
of the NDC in tunnel diodes. The DC voltage has to be 
added to the normal phase above Tc and the simplest ex- 
periment which we suggest is to measure the temperature 
dependence of the AC current (differential conductivity) 
as a function of the temperature. The theory reliably pre- 
dicts annulation of the differential conductivity at cooling 
below Tc- What exactly will happen at further cooling is 
difficult to be predicted experimentally: it could be NDC 
if the space homogeneity is conserved, creation of layered 
domain structure of normal and superconducting layers 
or complicated dynamic phenomena. Different instabil- 
ities can be triggered by small perturbations including 
the nature of contacts and chemical treatment of the sur- 
faces. One thing is sure: indispensably, new physics will 



happen in the area below the annulation of the differen- 
tial conductivity. The realization of the simplest possible 
scenario of creation of NDC will open the technological 
perspective to create a new type of high-frequency oscil- 
lators operating in the THz gapJ^ While in Ref. 13 is 
analyzed the case of bulk superconductors in the present 
work we demonstrated that for electric currents applied 
in the "dielectric" c-direction the dissipated power could 
be orders of magnitude smaller and perhaps Josephson 
coupled superconductors can realize the first technical 
applications of NDC in superconductors; it was one of 
the motivations to perform the present work. 

Few words we wish to add concerning the theoretical 
notions used in the present work. It seems very strange 
that Boltzmann kinetic approach has still limited usage 
in the physics of paraconductivity as it was for the normal 
conductivity hundred years ago. Just Green functions 
of the Boltzmann equation give the transparent picture 
what is happening in homogeneous electric fields. Sepa- 
ration of the variables in optical gauge with zero electric 
potential f = reduces the kinetic equation to an ordi- 
nary differential equation. Returning to the kinetic mo- 
mentum space gives the drift of the fluctuation Cooper 
pairs created by the electric field. 

Let us continue the discussion of physics in supercool- 
ing regime. Near to the annulation point of the differen- 
tial conductivity the paraconductivity is not a small ex- 
cess perturbation but will become comparable with the 
fluctuation conductivity. In the suggested experiment of 
applied DC electric field and small AC voltage one can 
observe in principle the AC component of the fluctua- 
tion magnetization M^. This physical effect is described 
by magnetoelectric susceptibility xme — dM^/dEz] this 
theory will be the subject of a future research. Measure- 
ments of magnetization and paraconductivity in one and 
the same sample can lead to exact determination of the 
time constant tq. This is the most important parame- 
ter of TDGL theory which reveals in part the nature of 
superconductivity and pairing mechanism. The 7r/8 mul- 
tipfier in Eq. ^ is only weak couphng BCS result, but 
as we already mentioned, for references of relevant works 
see Sec 4.2 of the review )^ the experimental accuracy now 
is not enough to observe reliably any deviations from this 
BCS 7r/8 value. There is no consensus that BCS theory 
is directly applicable for high-Tc cuprates and one might 
expect some correction multiplier Tj-ci to the life time of 
the order of one. The relative life time can give important 
information for the pairing mechanisms in cuprates. 

From qualitatively point the numerical value of tq is 
not so important because only fixes the scale for the 
graphical presentation of the experimental data. In the 
chosen dimensionless units the current voltage curves 
should be universal and the most important property 
of the current voltage curves is the predicted annula- 
tion of the differential conductivity at supercooling below 
Tc- Demonstration of this annulation and possible NDC 
in out-of-plane transport can trigger significant technical 
applications^'^ and we hope that experimental search of 



NDC in superconductors can start in the nearest future. 
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